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Abstract 

We consider the motion of a test particle in a one-dimensional system of 
equal-mass point particles. The test particle plays the role of a microscopic 
"piston" that separates two hard-point gases with different concentrations 
and arbitrary initial velocity distributions. In the homogeneous case when 
the gases on either side of the piston are in the same macroscopic state, we 
compute and analyze the stationary velocity autocorrelation function C(t). 
Explicit expressions are obtained for certain typical velocity distributions, 
serving to elucidate in particular the asymptotic behavior of C(t). It is shown 
that the occurrence of a non- vanishing probability mass at zero velocity is nec- 
essary for the occurrence of a long-time tail in C(t). The conditions under 
which this is a t~ 3 tail are determined. Turning to the inhomogeneous system 
with different macroscopic states on either side of the piston, we determine 
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its effective diffusion coefficient from the asymptotic behavior of the variance 
of its position, as well as the leading behavior of the other moments about 
the mean. Finally, we present an interpretation of the effective noise arising 
from the dynamics of the two gases, and thence that of the stochastic process 
to which the position of any particle in the system reduces in the thermody- 
namic limit. 
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I. INTRODUCTION 



Exactly solvable "toy" models are important in understanding the dynamic behavior of 
complex systems made up of a large number of particles. They also allow us to establish 
and understand the limitations of the approximations used in general to deal with systems 
of interacting particles, such as the Boltzmann equation. One such model, consisting of 
identical hard-point particles moving on a line and interacting through elastic collisions, 
was introduced several decades ago JO. Based on the observation that the particles merely 
exchange velocities in a collision, Jepsen 0] was able to calculate explicitly several properties 
of a gas of such particles. Subsequently, Lebowitz and co-workers |§- || refined and extended 
these calculations to include, among other aspects, a comparison with the results of the 
Boltzmann approximation. Very recently, this model system has been revisited by Piasecki 
and Gruber ||, 0, their main motivation being the construction of a one-dimensional analog 
of the "adiabatic piston" ||, with a central particle playing the role of a piston separating 
gases at different temperatures and densities to its left and right. Attention has also been 
focused on several related models from the point of view of the applicability of the Fourier 
law for heat flux in one- dimensional systems ||. Although the main interest in these specific 
contexts is in the case of an arbitrary ratio of the masses of the piston and a gas particle, the 
equal-mass case being a singular one in some sense, the analytical tractability of the latter 
makes it a valuable "theoretical laboratory" from which much can be learnt. In particular, 
it offers a model in which (i) one can pass to the thermodynamic limit in a rigorous manner, 
and (ii) the effect of recollisions can be completely and exactly taken into account. 

The renewed interest in the model prompts us to revisit it and extend its analysis. We 
focus on the motion of the central particle (which we shall refer to as the "piston" for 
brevity). A brief recapitulation of the relevant results from earlier work is given in Section 
II. In Section III we compute and study the stationary velocity autocorrelation function 
C(t) of the piston in the homogeneous case, when the gases on either side of the piston 
are in the same macroscopic state (i.e., equal densities and arbitrary but identical velocity 
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distributions), to amend and extend earlier results. The explicit expressions obtained for 
certain archetypical velocity distributions help us analyze the asymptotic behavior of C(t) 
to determine exactly when a power-law decay may be expected, and when the latter is a 
t~ 3 tail. Turning to the inhomogeneous system with different macroscopic states on either 
side of the piston, in Section IV we determine its effective diffusion coefficient from the 
asymptotic behavior of the variance of its position, as well as the leading behavior of the 
other moments about the mean. Finally, in Section V we show that there is an appealing 
and direct interpretation of the effective noise to which the many-particle interactions reduce 
in the thermodynamic limit as far as one-particle dynamics is concerned, and thus of the 
stochastic process represented by the position of any particle in the system. 



II. RECAPITULATION OF EARLIER WORK 

For ready reference, we record briefly the relevant results of earlier work that are needed 
for what follows, using the convenient notation employed by Piasecki |7j . One starts with a 
system of a finite number of particles, located in the interval [-L, L] of the x-axis, with free 
boundary conditions. The piston is initially located (without loss of generality) at X = 
and has an arbitrary initial velocity Vq. The iV~ particles to its left and N + particles to its 
right are initially at uniformly distributed random positions Xj (j = —N~, ... , 0, . . . , N + ), 
with independent, identically distributed velocities Vj drawn from normalized velocity dis- 
tributions ^(V) respectively. As the particles merely exchange velocities upon collision, at 
any instant of time the piston is on one of the "free" trajectories Xj+Vjt. This, together with 
the fact that the particles cannot move across each other, suffices to solve the problem of 
determining, among other quantities, the phase space distribution function of the piston (or 
that of any of the other particles) at any time t > by averaging over the initial positions 
and velocities of the gas particles on both sides of the piston ||, ||. In the thermodynamic 
limit L —>■ oo and — > oo such that one has finite densities lim N^/L = n on the left 
and right of the piston, the conditional one-particle distribution function of the piston is 
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found to be 

P(X,V,t\0,V o , 0) = 

e -t[n-a(v )+n+p(v )] Io (2t[n~ a(V )n + P(V )} 1/2 ) 5{X - V t) 5(V - V ) 
+ e -t[n-a(x/t)+n+p(x/t)] j [ n - ^ (y)6{Vt - X)9(X - V t) 
+ n + (/) + (V)6(X - Vt)0(V o t - X)] J {2t[n- a (X/t)n + (3(X/t)} 1/2 ) 



+ 



n-<f)-{V)6(Vt - X)0{V Q t - X) (n + [3{X/t)/n~a{X/t)) 1/2 + n + <f) + {V)6{X - Vt) 
x6{X-V t) (n-a{X/t)/n + p{X/t)) 1/2 ] I x (2t[n~a(X/t)n + /3{X/t)] 1/2 )} , 



(1) 

where 

oo W 

a(W) = J dU<p~(U)(U -W) , P{W)= J dU(f) + {U){W -U) , (2) 

W -oo 

and I n is the modified Bessel function of order n. For Vo = 0, this is the result recorded in 
Ref. 0. 

The velocity distribution of the piston is obtained by integrating P(X, V,t\0, Vo, 0) over 
X, and is given by 

P(V,t\0,V o , 0) = e -*[»-«W-^Wo)]j Q (2t[n- a (V )n + (3(V )} 1/2 ) 5(V - V ) 

/oo 
d jy e -t[»-<*(w)+n+/3(wo] | [ n -^-(y)0(y _ jy)0(jy - Vo) 
-oo 

+ - v)0(vb - wo] / (2t[n- a (iy)n + /5(iy)] 1/2 ) 

+ 

x 6(W - Vq) (n-a{W)/n + (3{W)) 1/2 ] h (2t[n~a{W)n + (3{W)} 1/2 ) J . (3) 

The stationary velocity distribution to which this tends as t — > oo is found by using the 
leading term in the asymptotic expansion of I n (z) for large argument, and carrying out the 
integral over the resulting Gaussian peaked at W = W, where W is the (unique) root of the 
equation 

n -a(W) = n + (3(W) . (4) 



n-<p-(V)d(V - W)9(V - W) (n + (3(W)/n~a(W)) 1/2 + n + 4> + (V)9(W - V) 



The normalized asymptotic velocity distribution of the piston is then found to be 

= n-<p-(Y)0(y ~ W) + n + <P + (V)6(W - V) 
E(W) ' 

where 

E(W) = -n-a'(W) + n + {3'(W) = n-\a'(W)\ + n + (3'(W) (6) 

is the normalization factor. In general, therefore, Poo{V) has a finite discontinuity at V — W. 
The asymptotic drift velocity of the piston, defined as (V(oo)) = dV V P^V) , is then 
trivially seen to be equal to W itself, on re-writing the definition of the latter in Eq. (Q) as 

oo W 

n~ J dUU(f)-(U) +n+ J dUU<p + {U) 
W=^ . (7) 

oo W 

nr J dU(p~(U)+n+ J dU<p+{U) 

W — oo 

In particular, for Maxwellian distributions (p^iV) at temperatures T ± , one has W ^ ac- 
cording as n~y/T~ ^ n + \ / T + . Since n(kBT ± ) 1 ^ 2 is essentially the linear friction coefficient of 
the corresponding gas, it is this dynamic property, rather than the pressure, that determines 
the direction of the drift of the central particle. The equality n~\jT~ = n + \/T + as the con- 
dition necessary for the piston to have a drift-free, purely diffusive motion asymptotically, 
has already been pointed out by Jepsen 0. 

We can proceed to show that, in the limit t —* oo, every particle attains the new station- 
ary velocity distribution Poo(V) (in the thermodynamic limit, of course). This remains true 
even if 0~ = <p + = <p, but n~ ^ n + . Thus, although the system under study does not have 
strong mixing properties, in the sense that the set of "free" trajectories {Xj + Vjt} is fixed for 
all time once the initial values {Xj} and {Vj} are specified, the effect of the collisions of the 
particles with the piston and among themselves is to mix the initial velocity distributions 
on the right and left of the piston. This happens by a "diffusion" of the "interaction front" 
of the piston with the other particles, through successive collisions. Only in the special case 
in which 

rT = n + = n, 0~ = + = , (8) 
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i.e., when the gases on both sides of the piston are initially in the same macroscopic state, 
does Poo(V) reduce to 4>(V) itself. We refer to this as the "homogeneous system" in what 
follows, in contrast to the more general inhomogeneous case. For simplicity, we shall also 
assume that <j>(V) = <j>(—V) throughout this paper. 

III. VELOCITY AUTOCORRELATION IN THE HOMOGENEOUS SYSTEM 

We now study the stationary velocity autocorrelation function of the piston, which re- 
quires that the initial distribution be a stationary one. As mentioned above, this only 
happens in the homogeneous case, to which we therefore restrict ourselves in the rest of this 
section. The results to be obtained are in fact valid for any particle in the system. 

Since W = in this case, the autocorrelation function is given by 

/oo poo 
dV / dV V VP(V, t\0,V , O)0O/ O ) • (9) 
-oo J — oo 

This quantity has been considered in earlier work both for a Maxwellian 4>(V) and in more 
general terms and been found to exhibit a t~ 3 power-law tail. However, the conclusions 
regarding the conditions under which this happens require modification; nor is the exponent 
of the power law invariably equal to —3. We therefore analyze the question afresh, extending 
and amending some of these earlier results. We also give a physical interpretation of the 
circumstances leading to a power law decay of C(t). 

On inserting Eq. (|3|) in Eq. ([J) and simplifying, we obtain the expression 

/>oo 

C{t)= / dW e~ nt[a(w)+ ^ w)] {/o (2nt[a(W)f3(W)} l/2 ) W 2 4>{W) 



+ nt (a(W) - Wa\W)f [ fc^ff A (2nt[a(WMW)} 1/2 ) 
-2I (2nt[a(W)P(W)} 1 / 2 )]} , (10) 

where the functions a(W) and (3{W) have been defined in Eqs. (H). It is readily verified 
that in the special case of a dichotomic velocity distribution 

<P + {V) = r(V) = <P(V) = ±[6{V + c) + 5(V - c)] , (11) 



one obtains from Eq. fllPl) the exponential decay 

C(t) = c 2 e~ 2nct . (12) 

Similarly, the known result for a Maxwellian 4>(V) is also recovered. 

C(0) is of course f dW W 2 (j)(W), and C(t) initially decreases linearly with t, with a slope 
that works out to 

/oo 
dW {W 2 (j)(W)a(W) + [a(W) - Wa'(W)} 2 } . (13) 
-oo 

In particular, for a Maxwellian <p(V) one can evaluate the integral involved to obtain 
(dC/dt)t=o = —^n^ksT/nra) 1 / 2 . In the general case, C(t) is a non-monotonic function 
of t, that becomes negative beyond a certain point and eventually approaches zero from 
below as t — > oo. The long-time behavior of the velocity autocorrelation yields valuable 
information on the mixing properties and memory effects in the system. The extraction 
of this behavior from Eq. ( [370] ) is non-trivial. It is helpful to note that ot(W) is a non- 
increasing function of W, with a(W) ~ — W as W — > — oo and a(oo) = 0; whereas (3{W) 
is a non- decreasing function of W, with /?(— oo) = and ^(W 7 ) ~ as — > oo. Further, 
whenever <j>~{V) = <p + (V) = <j>(V), we have (3{W) = a(-W) = W + a(W). An adequate 
number of terms in the asymptotic expansions of the Bessel functions and the other terms 
in the integrand in Eq. fllTf ) must be retained, consistent with the fact that non-vanishing 
contributions to the integral come from the region W 2 t < 0(1). As already mentioned, it 
has been shown 0, || that C(t) has a leading asymptotic behavior ~ t~ 3 for a Maxwellian 
<p(V) jllj. We have corroborated this, and also extended the result to the next term in the 
asymptotic expansion: after a very lengthy calculation, we find 

1 / m \ ( _ 5 \ 1 / m \ ( 3157T 367 



(ntf\2Txk B T) \} 2ir) 8(nt) 4 \2irk B T ) \ 177 16 tt 

(14) 

Based on the emergence of a t -3 tail in the Maxwellian case (and also for another extended 
distribution <p(V) that falls off like |l^|~ 3 ), as opposed to an exponential decay (Eq. (|12D ) for 
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the dichotomic distribution, it has been concluded |3j that C(t) decays exponentially if 0(V) 
is a compact ("finite") distribution, and has a power-law tail whenever the support of (f>(V) 
is non-compact. However, as we now proceed to show, the actual criterion for the emergence 
of a power law tail in C(t) turns out to be the existence of a non-zero probability mass at 

V = in 4>{y). (The physical reason for this will be described subsequently.) The issue 
is most clearly elucidated with the help of two simple yet archetypical distributions 4>(V) 
for which C(t) can be determined analytically, and exact asymptotic expansions obtained. 
These are, respectively, (i) a uniform distribution in the finite interval — c < V < +c, and 
(ii) a discrete distribution consisting of 5-functions at V — ±c and an additional one at 

V = 0. The former has a compact support, and yet C(t) turns out to have a t~ 3 tail; while 
the latter, in contrast to the dichotomic distribution of Eq. ([11]), leads in fact to an even 
heavier tail (~ t~ 3 / 2 ) for any non-vanishing weight of the central 5-function. We also extend 
the result in Eq. (|14"D to the case of a general 4>(V) that is sufficiently smooth at the origin. 

(i) Accordingly, let us consider the uniform distribution 



For this distribution, a(W) is respectively equal to —W for W < — c, and (c — W) 2 /4c for 
-c < W < +c; it vanishes identically for W > c. Recall also that (3(W) = W + a(W). 



closed- form expression for C(t). In terms of the dimensionless time r = net, this reads 



where \Fi is the usual confluent hypergeometric function. The slope at the origin is 
(dC(r)/dr) T= o = ~\° 2 - The analytic form of C(r) enables us to write down its exact 
asymptotic expansion for large r: 



4>(v) = ^-e(v + c )e(c-v) . 



(15) 



Inserting these in Eq. ( |T0D and carrying out the necessary calculations, we finally obtain a 




(16) 




(17) 



n=3 
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As mentioned earlier, this starts with an 0(t~ 3 ) term, although the support of 4>(V) is 
compact. Figure [V] depicts the long-time behavior of the correlation function, 
(ii) Next, consider the discrete distribution 

(f>(V) = 11 S(V) + — ^ [6(V + c)+S(V-c)), < n < 1 (18) 

which is an extension of the dichotomic distribution of Eq. ([TT| ) to include an additional 5- 
function at V = with a weight fi. Once again, c^VT) vanishes identically for W > c. For 
W < c , it is piecewise linear, being given by 



a(W) 



-W, W<-c 

|(l-/i)c-i(l + //)W, -c<W<0 (19) 
\(1-H)(c-W), 0<W<c. 

We obtain in this case (with r = net as before) 



C(r) = c 2 (l-n)e- T I l + (l-/i)r ^ dMe"' 1 "^ 



- — yi^ii (r<7(u)) - Jo (r</(u)) 



(20) 

1 /2 

where ^(tt) = ((1 — [i)(l — u)[(l — fi) + (1 + jJ>)u]) . The slope at the origin is 
(dC(r) I 'dr) T=0 = —(1 — /i)(2 — /i)c 2 . The long-time behavior in this case is however quite 
different from that found in the previous cases. Owing to the singularity in <p(V) at the 
origin, a'(W) has a jump at W = 0. As a consequence, the integrand in Eq. ( fL0|) is now a 
function of \W\ rather than W. This leads to the occurrence of both even and odd powers 
of | W | in the small- W expansion of the integrand, using which the asymptotic expansion of 
C(t) is determined. For the latter, we now obtain 

rM /ic 2 (1-Zi) 1 / 2 /ic 2 3(3 -2/i 2 ) 
{) ~ r 3 / 2 4 (2tt) 1 /2 r 5 / 2 32 (2 tt) 1 ^ (1-^)1/2 • ^ > 

Thus C(t) now has an even slower power-law decay, starting with an 0(t~ 3 ^ 2 ) term, as long 
as n 0, i.e., as long as there is a finite probability mass at V = 0. When fi = 0, all the 
terms in the asymptotic expansion vanish, and C(t) reverts to the exponential decay that 

10 



obtains in the case of the dichotomic distribution, Eq. ([12]). Figure [V] shows the long-time 
behavior of the correlation function for different values of the weight parameter /x, including 
(for ready comparison) the case \x = 0. 

The relative roles of the 5-functions in <f>(V) at V = and at V — ±c may be examined 
a little more closely. This aspect is not so transparent in the representation of Eq. fl2"DD for 
C(t), but is made more manifest with the help of its Laplace transform. This enables us to 
write C(t) in the form 



C(t) = (l-fi)c 



2 



_ 2net 

e !+f 

1 + n 



, , r -i (sy+[s* + 2snc(l 1 

^ ((l+fi)s + 2nc)(s + fx[s 2 + 2snc(l-^)} 1 / 2 + 2nc(l- fi)) j ' 1 ' 

where £ _1 denotes the inverse Laplace transform. Comparing this with the pure exponential 
decay c 2 e~ 2nct that obtains for the dichotomic velocity distribution, we see that the 5- 
function in 4>(V) at V = is entirely responsible for the second term (which vanishes 
when fi = 0). Further, the timescale in the exponential part is itself modified from the 
usual correlation time for a dichotomic process, which is (2nc) _1 in the present context, to 
(1 + n)(2nc)' 1 , as one might expect on physical grounds. 

We are now in a position to understand the physical origin of the power-law tail in C (t) . 
The particles of the system do not undergo any systematic drift in the homogeneous case. 
Going back to an inspection of the manner in which the particle under consideration skips 
from one free trajectory to another through collisions, we see that, if the stationary velocity 
distribution <p(V) of the gas particles has a finite probability mass at V — 0, the particle will 
repeatedly find itself on a trajectory with zero slope, i. e., revert to the zero (= average) 
velocity state. This persistence is like a memory effect, and it shows up as a slow (power-law) 
decay of C(t). The compactness or otherwise of the support of <j)(V) does not play a role as 
far as this aspect is concerned. 

It is also possible to find the precise conditions under which the leading asymptotic 
behavior of C(i) starts with a t~ 3 term: this is so if <f)(V) is at least twice differentiable at 
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the origin, and moreover 0(0) and (j)"(0) do not both happen to be zero. (We recall that 
4>(V) has been taken to be a symmetric function, so that all its derivatives of odd order 
vanish at the origin.) The general asymptotic expansion of C(t), for a distribution (f>(V) 
which is differentiate a sufficient number of times at V = 0, reads: 

C(t) ~ [0(0) - 6a(O)0 2 (O) - « 2 (O)0"(O)] 

' [-315 0(0) + 3456 a(O)0 2 (O) - 2880 a 2 (O)0 3 (O) 



256 a(0) {ntf 

+ 19840 a 3 (O)0(O)0"(O) - 2208 a 2 (O)0"(O) - 256 a(O)0 ( ™ } (O)] - 



(23) 



This extends the result presented in Eq. (|l4j) for a Maxwellian. Thus, for a distribution 
4>(V) that is regular at the origin and has a non- vanishing derivative of some finite order at 
that point, implying that there is a non-zero probability mass at V — 0, C(t) will certainly 
have a power-law decay: If 0(0) ^ 0, the leading term is generically ~ t~ 3 ; on the other 
hand, if 0(0) = and its first non- vanishing derivative at the origin is its (2r)-th derivative, 
the leading term in C(t) is ~ t~ r ~ 2 . 



IV. THE INHOMOGENEOUS SYSTEM 

We turn now to the inhomogeneous system, in which the particles to the left and right 
of the piston are initially in different macroscopic states specified by (n~,(f>~) and (n + ,0 + ) 
respectively. The piston now has, in general, a non-vanishing mean drift velocity that 
asymptotically approaches W. However, as we shall see, the variance of its position indeed 
increases linearly with time. The quantity of interest is therefore the effective diffusion 
coefficient D, which we shall determine. We also find the leading asymptotic behavior of 
the other moments of the position (about its mean value). 

In the homogeneous system, D is of course equal to §^C{t)dt (this integral being 
absolutely convergent for the system at hand). However, owing to the non-stationarity of 
the velocity autocorrelation in the inhomogeneous case, D must now be computed directly 
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from the long-time behavior of the mean square displacement of the piston. The asymptotic 
behavior of the piston does not depend on its initial state. We can therefore set Vq = in 
P(X, V, t\0, Vo, 0) (Eq. ([[])) and integrate it over V to calculate the position distribution 
function p(X, t) of the piston. Using the fact that the derivatives of the functions a and (3 
are given by 

oo W 

a'(W) = - J dU<f)-{U) and (3'{W)= JdU0 + (U), (24) 

W -co 

we find 

p(X,t) = e -*h-«(o)+" + /?(°)]j (2t[n~a(0)n + {3(0)} 1/2 ) 5(X) 

+ e -t[n-a(X/t)+n+0(X/t)] j [ n -Q(X) \a\X/t) | + n + 6 (- X) ft (X / 1)] 

x I (2t[n~a(X/t)n + f3(X/t)f 2 ) + \rTB{-X)\a'{X/t)\ (n+ (3{X/t) jn a[X /t)f' 2 

1/2 - 



+ n + 9{X)(3'{X/t) (n~a{X/t)/n + (3{X/t)y' 2 h (2t[rT a{X/t)n + (3{X /t)} 1 / 2 ) 

(25) 

The variance of the position is given by 

/oo 
dX (X - (X(t))f p(X,t) . (26) 
-oo 

In the long-time limit, (X(t)) = Wt. Using the asymptotic behavior of the Bessel functions 
in p(X, t), the leading behavior of the variance is given by 

4-5/2 roo 

((X-Wt) 2 )~- -=-—[n-\a'{W)\+n + (3'{W)] / dW(W - W) 2 

xexp\-t(W -W) 2 ([n-/Aa(W)} 1/2 \a'(W)\ + [n + /A(3(W)} 1/2 (3'(W)) 2 } (27) 



which simplifies to 2Dt, with a diffusion coefficient given by 

rra{W) 1 n-a{W)+n + (3(W) , . 

J-* — ^ ^ 9 — ^ ^ 9 • \^°) 

[n-\a'(W) | + n+(3'{W)] 2 [n-\a'{W) | + n+P'{W)] 
This is the general formula sought. 
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We first note (as a check) that in the special case of the homogeneous system, Eq. ( [28] 
reduces to 



oo 



D= ^i = if duum = m, (29) 

n n } 2n 



o 



in agreement with the known result ||. As mentioned earlier, in this case D must also 
be equal to the integral of the velocity autocorrelation C(t). We have verified that this is 
indeed so. 

Some interesting special cases emerge from the general formula of Eq. (p8|). If the 
densities n~ and n + are such that the drift velocity W = even though one has different 
(but symmetric) distributions <fi + {V) and 4>~{V) on either side of the piston, the formula for 
D simplifies somewhat. Since o/(0) = — /3'(0) = —1/2, we find 



n «(0) 
(n~ + n + y 



D = —- I iL ■ (30) 



In particular, if ± are Maxwellians (with n \[T~ = n + \/T + to ensure that W = 0), 

n- {8k B T-^ 1/2 



D = ——— (31) 
(n~ + n + ) 2 \ Tim J 

On the other hand, if <j)~{V) = <f) + {V) = <j){V), but the system remains inhomogeneous 
because n + ^ n~, then W ^ 0. For the compact uniform velocity distribution of Eq. (|15|). 
we find 



— = c v n vn (32) 
'n + v n H 



and 



D = = = . (33) 
(\/n~ + yn + ) 2 



For the discrete distribution of Eq. ( 18l) , we have 



TI7 (l-/x)(n -71+) 

w = c — t; \ 7^ c > (34) 

n> 1 - n + n< 1 + fi 
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where n> = max(n , n + ) , n < = min(n , n + ). The corresponding diffusion coefficient is 
found to be 

n = 4cw-ra+(l -//) 

Mi-zil + Mi + z*)]'" 

It is noteworthy that the interplay of the central 5-function in the velocity distribution and 
the inhomogeneity due to the different densities on either side of the piston affects even the 
diffusion coefficient. Setting \i = in the above yields the corresponding expressions for the 
dichotomic distribution of Eq. ( |iTl) . 

The leading asymptotic behavior of the higher moments ((X — W t) r ) can also be de- 
termined. Here we merely quote the salient result obtained. For the even moments r = 21, 
it is straightforward to show (along the same lines as in the case of the variance) that 
((X — Wt) 21 ) ~ 0(t ). The calculation is more involved for the odd moments r = 21 + 1, 
but the final result is that {(X — W t) 2l+1 ) ~ 0{t l ) as well. As the expressions obtained for 
the precise coefficients are lengthy, we do not write them down here. 



V. INTERPRETATION AS A STOCHASTIC PROCESS 

We conclude by showing that (in the thermodynamic limit) the form of the distribution of 
the position of the piston, in fact that of any of the particles in the system, is effectively that 
of a stochastic process driven by a noise that can be given a direct physical interpretation. 

Conventionally, the stochastic approach to single-particle dynamics in a many-body sys- 
tem begins with its modeling by a stochastic evolution equation involving noise terms with 
prescribed statistical properties. One then extracts the corresponding properties of the 
driven variable(s). Here, however, we have the converse situation. The exact time- dependent 
one-particle distributions are known, and the task is to identify the stochastic process to 
which the complicated dynamics effectively reduces, at least as far as the one-particle dy- 
namics is concerned. What kind of stochastic process does the position X(t) of the piston (or 
any other particle) represent, after the averaging over the initial states of the gas particles is 
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done, and the thermodynamic limit taken? And in what kind of "noise" are the combined 
effects of the other particles in the system encapsulated? 

It is evident from the rather complicated expressions for P(X,V,t\Q,Vo, 0) and the re- 
duced distributions derived from it that X(t) is unlikely to satisfy any simple or standard 
stochastic differential equation; nor does p(X, t) appear to be the solution of any simple 
master equation - in particular, of any obvious partial differential equation of finite order. 
Intricate correlations exist, that cannot be neglected. The effects of recollisions are ob- 
viously significant, a direct instance being provided by the form of the first term in Eq. 
(HI). This term represents the probability for the piston to find itself in its initial state at 
time t. Now, the probability that the initial state of the piston (X = X = 0, V = Vq) 
persists till time t (i.e., the piston suffers no collisions till time t) is easily shown to be 
simply exp [— t(n~a(Vo) + n + (3(Vo))}. Thus the extra factor I in the term proportional to 
5(X - V t) 8{V - Vq) in P(X, V,t\0, V , 0) is entirely due to the effects of recollisions 0. 
As the concept of an effective noise is only meaningful in the thermodynamic limit and 
when ergodicity obtains, we must examine for this purpose the structure of the terms in the 
solutions other than the ones arising from the returns to any specific initial state. 

The occurrence of the Bessel functions J and Ji in Eqs. (|I|) and Q25| ) seems to suggest 
some sort of link with dichotomic diffusion (i.e., the integral of a dichotomic Markov process) 
and the well known telegrapher's equation and its solution. Indeed, in the homogeneous case, 
with 4>(V) equal to the dichotomic velocity distribution of Eq. (|TTD and V — ± c , Eq. fl25|) 
for p(X, t) does reduce to the solution correponding to dichotomic diffusion [pi]] , once again 



except for the extra factor of Io(nct) in the "ballistic" term representing the probability of 
the occurrence of the initial state at time t. But this does not explain the origin of the Bessel 
functions in the general case. Nor does it really do so even in the special case referred to, 
other than the not-very-helpful observation that the solution of the telegrapher's equation 
involves Jo and I\. As the effective "noise" we seek should be essentially the same for every 
particle, our arguments should indeed apply to any of the particles, and not just the piston. 
Proceeding as in the case of the piston (i.e., averaging over the initial positions and velocities 
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of all the particles except the piston, and with X = 0, Vq = 0), we find the following result 
for the position distribution of the b th particle at time t: 

p b (X,t) = e -i n ~ a{0)+n+m ] (n + p(0)/n-a(0)) b/2 I b (2t[n-a(0)n + f3(0)} 1/2 ) 5(X) 
+ e -t[n-a(x/t)+n+f3(x/t)] (n+ /t) / n~ a(X / t)) b/2 

x {(n-0(X)\a'(X/t)\ + n + 6{-X)(3'{X/t)) I b {2t[rCa{X/t)n + P(X/t)] 1 ' 2 ) 

+ rre(-X)\a'(X/t)\ (n + /3{X/t)/n-a{X/t)) 1/2 I b+1 (2t[ n -a{X/t)n + (3{X/t)] 1/2 ) 

+ n + 6(X)(3'(X/t)(n-a(X/t)/n + (3(X/t)) 1/2 J 6 _a {2t[n~ a(X/t)n + P(X/t)} 1/2 ) } . 

(36) 

Here b G Z. (Setting 6 = one recovers the result in Eq. ( plf) for the piston, remembering 
that J_i = The occurrence of the Bessel functions h+i, h and I b -i shows quite clearly 
that we must look for a link to Bessel functions other than one via the solution to the 
telegrapher's equation. 

The example of the dichotomic velocity distribution in the homogeneous case does pro- 
vide a valuable clue, though. Let us therefore examine this case for a moment, focusing 
on the piston, and taking Vo also to be either c or — c, to mask the effects of any special 
initial conditions. The actual "free" trajectories of all the particles are then straight lines 
with slopes restricted to the values ± c. A little thought shows that after its first collision, 
the piston alternately 'rides' on a free trajectory belonging to the gas on its right, and one 
belonging to the gas on its left. For brevity, we shall refer to these as 'right' and 'left' tra- 
jectories. (Which of these the piston gets on to first, depends on whether its initial velocity 
Vq is equal to c or — c.) Moreover, the piston is alternately hit by a 'right' particle with 
velocity — c and a 'left' particle with velocity +c. The number of 'right' collisions minus the 
number of 'left' collisions can only take on the values +1, and —1. The resulting zig-zag 
path is precisely that of a particle whose position X satifies the stochastic differential equa- 
tion X = c£(t), where is a stationary dichotomic Markov process (DMP) alternating 
between the values ±1 with a certain mean switching rate A. Such a DMP is generated by 
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a stationary Poisson pulse process of intensity A. One could also regard it as made up of 
two independent Poisson pulse processes, each with an intensity A/2, alternating with each 
other. This would seem to be a little more closely linked to the present situation, where 
one might imagine the two states of £(£) to be related in some sense to the piston being 
on a right trajectory and a left trajectory, respectively. But the connection is still far from 
obvious, and requires some more work. 

Let z/ + and v~ be two independent stationary Poisson processes with respective intensities 
(i.e., mean rates) A + and A - , so that their mean values at time t are A + 1 and A - t. It is 
easily shown that their difference (u + — u~), which can take on any integer value, has a 
time-dependent distribution given by 

Pr (v+ - v- = r; t) = e ~( A++A ^ (A+ t) r/2 I r hy/(X+t) (A"t)) , r E Z . (37) 

It is this distribution that holds the key to understanding the structure of the one-particle 
distributions in the problem under consideration. Let the piston be at a position X > 
at time t, on a (segment) of a free trajectory. Translating the entire system to bring the 
initial coordinate of this trajectory to the origin, the instantaneous velocity of the piston 
is X/t (recall that the trajctories are all straight lines, a direct consequence of the equal 
mass condition). It can be hit by a left particle provided the latter has a positive velocity 
(+c in the case under consideration) that is greater than X/t. The mean rate at which 
this happens is given by the product of the concentration n of the gas, the magnitude of 
the relative velocity (c — X/t), and the probability | that the velocity of the gas particle 
is c (see Eq. (|TT|) ): in other words, A - = \ n(c — X/t). Similarly, the mean rate at which 
the piston is hit by a right particle is given by A + = \n(c + X/t). Moreover, the number 
of 'right' collisions minus the number of 'left' collisions only takes on the values +1, and 
— 1. Putting in all the foregoing facts and their obvious extension to the case X < 0, and 
transforming from the random variable (z/ + — v~) to X, we are led to the expression 

|e-" rf {/„ (n(c 2 t 2 - X 2 ) 1/2 ) + {{ct + X)/(ct - X)) 1/2 h (n(c 2 t 2 - X 2 ) 1 / 2 ) } 

x 9(X + ct) 6(ct - X) . (38) 
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This is precisely the solution for p(X, t) to which Eq. (p5|) reduces in this special case, apart 
from the contribution from the initial state. The latter is e~ nct 5(X ± ct) when Vq = ± c, 
and e~ nct Io (net) S(X) when Vq = 0. It is also worth noting how the factors of 2 (coming 
from the formula of Eq. fl3~T|)) and | (coming from the rates A 1 * 1 ) cancel out in the argument 
of the Bessel functions in (|38|). 

These arguments are extended to the general inhomogeneous case as follows. When 
the piston is at position X at time t, collision with a left particle is possible provided the 
latter has a velocity U in the range (X/t, oo), and the magnitude of the relative velocity is 
(U — X/t). Since the gas on the left has a concentration n~ , and the velocities of its particles 
are drawn from the distribution <p~ , the effective mean rate of 'left' collisions of the piston 
is given by 

X~ — n~ J dU(jr{U){U - X/t) , (39) 
x/t 

which is nothing but n~a(X/t) . Similarly, in the same given state the piston can only be 
hit by a right particle with a velocity in the range (—oo, X/t), and it follows that 

x/t 

A+ = n + J dU<p + (U)(X/t - U) = n + f3(X/t) , (40) 

— oo 

since it is the magnitude of the relative velocity that appears in the mean collision rate. 
This explains the genesis and form of the argument of the Bessel functions in Eq. (^). The 
extra factors involving n~\a'(X/t) \ and n + /3'(X/t) that appear in the expression for p(X, t) 
in Eq. (|25f ) are just the Jacobians that arise when we transform from the distribution of 
(v + — v~) to that of X. Finally, although the piston can have successive left collisions or 
right collisions in the general case, unlike what happens in the case of the dichotomic velocity 
distribution, the remarkable fact is that their contribution to the probability distributions 
seems to vanish in the system at hand. The number of 'right' collisions minus the number of 
'left' collisions only takes on the values +1, and —1 even in the general case, presumably as 
a consequence of the smearing out implied by the averaging and the thermodynamic limit. 
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Likewise, for the 6 th particle from the piston, the difference between the number of times 
the particle has occupied a right trajectory and the number of times it has occupied a left 
trajectory up to any time t predominantly takes on the values b and b ± 1 in this limit. The 
central value b of this difference is in fact a consequence of its initial ordering, that is not 
altered by the collisions. The foregoing provides an understanding of the structure of the 
terms (other than the contribution due to the initial state) in Eq. (|36|). 

In conclusion, we see that the motion of any particle in the system may be regarded, 
in the thermodynamic limit, as being driven by two independent Poisson pulse processes 
representing the effects of the gases on the left and right of the central particle. The in- 
tensities (mean rates) of these processes have the direct physical interpretation given above 
with regard to Eqs. (|39D and (|40|). As a Poisson process per se is an uncorrelated pulse 
process, each particle is effectively subjected to two independent noises in this precise sense. 
However, as the intensity of each noise is state-dependent (the driven variable X appears 
explicitly in the limits of integration in a(X/t) and (3(X/t)), the flow of X is not given by 
any simple stochastic differential equation with additive or even multiplicative noise, which 
is only to be expected. 
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FIGURES 

FIG. 1. Asymptotic behavior of the normalized velocity autocorrelation C(t)/C(0) as a func- 
tion of time (in units of 1/nc), for a uniform distribution (j)(V) (Eq. 



FIG. 2. Asymptotic behavior of the normalized velocity autocorrelation C(t)/C(0) as a func- 
tion of time (in units of 1/nc) for different values of //, for the discrete distribution (j)(V) in Eq. 
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